Theoretical issues with AS in three-body B decays 
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Theoretical uncertainties with the time-dependent CP asymmetries AS and Acp in B c — > 
K+K-Ks, K S K S K S and K s tt°tt decays are discussed. In order to have a reliable estimate of 
CP asymmetries, it is very crucial to understand the underlying mechanism for the nonresonant 
background which dominates the KKK modes. 
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I. INTRODUCTION 

Time-dependent mixing-induced and direct CP asym- 
metries Sf and Af, respectively, in charmless three-body 
decays of the B mesons have been measured by BaBar 
and Belle for K S K S K S and R + R~R° final states. In 
addition, BaBar has also measured time-dependent CP 
asymmetries for K + K~ Kl and Rstt q tt° and performed 
a Dalitz plot analysis for K + K~ R° . 

For CP eigenstates, the measured mixing- induced CP 
asymmetry Sf can be used to define an effective sin 2/3 
via 



s f = ~Vf sin 2 A=ff 



(i) 



with 



rjf being the CP eigenvalue of the final state /. 

" o 



While K S K S K S , K S K S K L and K S 7r°7r u have fixed CP- 
parities with rjf = 1, —1, 1 respectively, R + R~R$ is an 
admixture of CP-even and CP-odd components, ren- 
dering its CP analysis more complicated. By excluding 
the major CP-odd contribution from 4>Rs, the 3-body 
K + K~Ks final state is primarily CP-even. A measure- 
ment of the CP-even fraction / + in the B° — > R + R~Rg 
decay yields /+ = 0.89 ± 0.08 ± 0.06 by BaBar and 
0.93±0.09±0.05 by Belle 0, while the CP-odd fraction in 
K+K-K L is measured to be /_ = 0.92±0.33+°^±0.10 
by BaBar 0. Hence, while rjf = 1 for the KsK$Ks 
mode, rj f = 2/+ - 1 for R + K~K S and r) f = -(2/_ - 1) 
for R+R-R L . 

The world average of sin 2/3 measured in the de- 
cays B — > 4>ccKs with (j> cS being a charmonium 
is sin2/3 0ce K s = 0.675 ± 0.026 0. However, the 
time-dependent CP asymmetries in the b — > sqq 
penguin-induced two-body decays such as B° — » 
(</>, w, 7T°, 77', fo)Ks measured by BaBar and Belle indi- 
cate sin2/3 ff = 0.53 ± 0.05 0. This seems to im- 
ply a 2.6 a deviation from the expectation of the stan- 
dard model (SM) where CP asymmetry in all above- 
mentioned modes should be equal to sin 2/30 ce ^ s with 
a small deviation at most 0(0. 1) 0. 

The measured results of sin2/3 c ff for R + R~ Rs,l, 
RsKgRs and Rstt°tt° are listed in Table IV below. In 
order to see if the current measurements of the deviation 
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of sin2/3 e ff in these modes from sin2/30 cE if s signal New 
Physics in b — > s penguin-induced modes, it is of great 
importance to examine and estimate how much of the 
deviation of sin2/3 ff is allowed in the SM. One of the 
major uncertainties in the dynamic calculations lies in 
the hadronic matrix elements which are nonperturbative 
in nature. One way to circumvent this difficulty is to im- 
pose SU(3) flavor symmetry or isospin and U-spin 
symmetries to constrain the relevant hadronic matrix 
elements. While this approach is model independent in 
the symmetry limit, deviations from that limit can only 
be computed in a model dependent fashion. In addition, 
it may have some weakness as discussed in 0. 

As mentioned above, BaBar has performed a Dalitz 
plot analysis for B° —> R + R~R° decays and found 
that it is dominated by 5-wave nonresonant contributions 
with fraction 91 ± 19% 0. A large nonresonant compo- 
nent is also found in B + — » R + R + R~ decays [ToL lll|. 
This is a surprise because it is known that in three-body 
decays of D mesons, the nonresonant background is at 
most 10%. Hence, it is very important to take into ac- 
count the large nonresonant effects in any realistic model 
calculations. 



II. FORMULISM FOR CHARMLESS 3-BODY B 
DECAYS 

We consider the decay P° — > R + R~R° as an illus- 
tration. Under the factorization approach, its decay am- 
plitude consists of three distinct factorizable terms: (i) 

the meson emission process, (B — > R + R°) x (0 — > R~), 

(ii) the transition process, (B — > R ) x (0 — > R + R~), 



and (iii) the annihilation process (B — > 0} x (0 — > 

R + R~R ), where (A — > B) denotes a A — > B transi- 
tion matrix element. 

For the kaon emission process, the amplitude induced 
by V — A currents reads 

A:urrent-ind = (K~ (p 3 ) | (su) V -A |0) 

x (R (p 1 )R+(p 2 )\(ub)v^A\B a ) 

JK rr) 2 , 1 2 2 \ 

= — —[2m 3 r + (m B - S12 ~ m 3 )ui + 
+ (s 2 3 - si 3 - ml + m\)u)-}, (2) 
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FIG. 1: Point-like and pole diagrams responsible for the B — > 
K + K matrix element of the current U7 M (1 — 75)6, where the 
symbol • denotes an insertion of the current. 



where = {pi + Pj) 2 , (qiq2)v-A = qil^ - 75)92, and 
r, cj+ and cj_ are three unknown form factors. In prin- 
ciple, one can apply heavy meson chiral perturbation 
theory (HMChPT) to evaluate the form factors r, ui + 
and uj- [121 ] . The relevant diagrams are depicted in 
Fig. 1. However, this will lead to too large nonreso- 
nant decay rates in disagreement with experiment [l3|. 
A direct calculation indicates that the branching ratio of 
B° — > K + K~K° arising from the current-induced pro- 
cess alone is already at the level of 77 x 10~ 6 which sub- 
stantially exceeds the measured total branching ratio of 
25 x 1CP 6 jij. The issue has to do with the applicabil- 
ity of HMChPT. In order to apply this approach, two of 
the final-state pseudoscalars have to be soft. The mo- 
mentum of the soft pseudoscalar should be smaller than 
the chiral symmetry breaking scale of order 1 GeV. For 
3-body charmless B decays, the available phase space 
where chiral perturbation theory is applicable is only a 
small fraction of the whole Dalitz plot. Therefore, it is 
not justified to apply chiral and heavy quark symmetries 
to a certain kinematic region and then generalize it to 
the region beyond its validity. 

To circumvent the aforementioned difficulty, we as- 
sume the momentum dependence of the nonresonant am- 
plitude in the exponential form 



.4 



current — ind 



^HMChPT 
current— ind ^ 



-a NR PB ■ (Pi +P2 ) p i<Pi2 



(3) 



so that the HMChPT results are recovered in the soft me- 
son limit pi,p2 — ► 0. That is, the nonresonant amplitude 
in the soft meson region is described by HMChPT, but 
its energy dependence beyond the chiral limit is governed 
by the exponential term e~ a ^ PB '^ Pl+P2 K The parame- 
ter a NR can be constrained by the decay B~ — ► 7r + 7r~7r~ 

M- 

For the kaon emission process, there exist two more 
amplitudes, one induced by currents and the other in- 
duced by scalar densities: 

M = (I?(p3)\(sd) V - A \0)(K+(p 2 )K-(p 3 )\(db) v _ A \B*), 



A 2 



-ftO, 



(lT(p 3 )\rsd\0){K + (p 2 )K-(p 3 )\db\B 



(4) 



OZI suppressed, they do receive intermediate vector me- 
son and scalar pole contributions, respectively. For ex- 
ample, 

(K+(p 2 )K-(p 3 )\(db) v _ A \B°) R = 



E 



9 



ir fo . - s 23 - vm fo . 



(f 0i \(db) v ^ A \B°), (5) 



where denote the generic /o-type scalar mesons, f Qi — 
/ (980), / (1370), / (1500), X (1550), ■ • • . 

For the transition amplitudes, we need to con- 
sider the two-body matrix elements {T^+i^ - 1 |0) and 
(K + K~\ss\0). Both receive resonant and nonresonant 
contributions. The first matrix element can be related 
to the kaon electromagnetic (e.m.) form factors F^ n K 
and K for the charged and neutral kaons, respec- 
tively. The nonresonant form factor can be constrained 
by the asymptotic behavior implied by pQCD [l5j |. The 
second matrix element has the expression 



(K+(p 2 )K-(p 3 )\ss\0) = 



m foJf 0i 9 



f 0i -*K+K- 



7o, 



S23 



l /o, r /o 



+ /. 



NB 



(6) 

is defined by 



Although the 3-body matrix elements in A\ and A 2 are 



where the scalar decay constant /| o 
(/oi|ss|0) = m /o ./; Q .. 

It turns out that the nonresonant background in B — » 
P\P 2 transition described by Eq. ([3]) it too small to 

account for the experimental observation that B — > 
K + K~K is dominated by the nonresonant contribu- 
tions. This implies that the two-body matrix element 
e.g. {KK\ss\0) induced by the scalar density should 
have a large nonresonant component. In the absence of 
first-principles calculation, we will use the decay mode 

B° — K S K S K S or B~ -> R-K S K S to fix f™ pi]. 

The calculated branching ratios of resonant and non- 
resonant contributions to B — > K + K~K and Kstt tt° 
are summarized in Tables Q] and UH respectively. The 
theoretical errors shown there are from the uncertainties 
in (i) the parameter a NR which governs the momentum 
dependence of the nonresonant amplitude [cf. Eq. (|3"|)]. 
(ii) the strange quark mass m s , the form factor F^ K 
and the nonresonant contribution f^ R constrained by 
the KsKsKs rate, and (iii) the unitarity angle 7. We 
see that the predicted rates for resonant and nonresonant 
components are consistent with experiment. The nonres- 
onant background arises dominantly from the transition 
process (90%) via the scalar density induced vacuum to 
KK transition, namely, (K + K~\ss\0), and slightly from 
the current-induced process (3%). In the Ksir°n° mode, 
the nonresonant fraction is found to be of order 40%. 
Recall that large nonresonant contributions to the de- 
cays B~ — » K~ir + Tr~ and B — > K ir + ir~, at the level 
of (35 — 40)%, have been reported by Belle [l6j |. 

The K + K~Kg mode is an admixture of CP-even and 
CP-odd components. By excluding the major CP-odd 
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FIG. 2: The K + K mass spectra for B° -> K+ K~ K s decay from (a) CP -even and (b) CP-odd contributions. The insert in 
(b) is for the <j> region. The full K + K~ Ks spectrum, which is the sum of CP-even and CP-odd parts, measured by BaBar 
is depicted in (c). 



TABLE I: Branching ratios (in units of 10 6 ) of resonant and 
nonresonant (NR) contributions to B — > K + K~ K . Theo- 
retical errors correspond to the uncertainties in (i) qnr, (h) 
m s , Fj? K and onr, and (iii) 7. Experimental results are taken 
from 



Decay mode BaBar [9] 



Theory 



4>-F 


3.07 


± 


0.45 


9 e+0. 0+0. 5+0.0 
^•0-0.0-0.4-0.0 


/o (980)^° 


5.31 


± 


2.19 


r o+0. 0+0. 1+0.0 
°' c5 -0. 0-0. 5-0.0 


A- (1550)A"- 


0.98 


± 


0.44 


n Qo+0. 00+0. 16+0. 00 
u - yo -0. 00-0. 15-0. 00 


NR 


21.7 


± 


5.1 


■\A /l+0. 4+2. 1+0.1 
14 - 4 -0. 5-1. 5-0.1 


total 


23.8 


± 


2.0 ± 1.6 


- i +u.3+y.ii+u.i 

1 1 - l -0. 3-2. 9-0.1 



TABLE II: Same as Table U except for the K°n°iT mode. 



TABLE III: Branching ratios for B° -> 
K+K'Ks, K S K S K S , K s KsK l Kstv°tt decays and the 
fraction of CP-even contribution to B -> K + K~K S . The 
branching ratio of CP-odd K + K~ Ks with <f>Ks excluded 
is shown in parentheses. Results for (K + K~ Kl)cp± are 
identical to those for (K + K~ Ks)cp^- Experimental results 
are taken from 



Final State 



K + K~K S 
{K+K-K s )cp+ 
(K+K-K s )cp- 

KsKsKs 
KsKsKl 



theory 



o r^+0. 13+1. 82+0. 06 

o-^^-0. 14- 1.45-0. 06 
c Q7+0. 04+1. 37+0. 04 
D.» I _o. 04-1. 12-0. 04 
1 C7+0. 09+0. 46+0. 02 
J-O ' -0.10-0.32-0.02 
(n 1/1+0.06+0.14+0.01 \ 
V u - ly: -0. 06-0. 06-0. 01 1 

input 

a na+0- 11+0.61+0.02 
°- uo -0. 08-0. 69-0. 02 
r, n/i+0. 01+1. 96+0. 03 
D - z *-0. 02- 1.45-0. 04 



g(10 



)cxpt 

12.4 ± 1.2 



■.2 ± 0.9 

< 14 



Decay mode 
/o(980)7? u 

1TV 

Kl°(U30)n° 
NR 



Theory 



total 



q o+0. 0+2. 0+0.0 
°-°-0. 0-0. 4-0.0 

n kk+o. oo+o. i6+o. oo 

u -°°-0. 00-0. 13-0. 00 
9 n+0. 0+0. 8+0.0 
z - o -0. 0-0. 6-0.0 
r -I +0.1 + 1.8 + 0.0 
"•- L -0. 1-1. 0-0.0 

+U.U+3. 9+0.1 



12.5 



-0.0-2.9-0.1 



contribution from (j>Ks, the 3-body K + K~K$ final state 
is primarily CP-even. The K + K~ mass spectra of the 
B° — > K + K~Ks decay from CP-even and CP-odd con- 
tributions are shown in Fig. [2j For the CP-even spec- 
trum, there are peaks at the threshold and tjik+k- =1.5 
GeV region. 1 The threshold enhancement arises from the 



fo(980)Ks and the nonresonant fg +K contributions. 
For the CP-odd spectrum, the peak on the lower end cor- 
responds to the 4>Ks contribution, which is also shown 
in the insert. The small middle hump in Fig.[^a) comes 
from the interference between b — > s and b — > u ampli- 
tudes. The b — ► u transition is governed by the current- 
induced process (B -> K+K } x (0 -> K~). From Eq. 
<j3j> it is clear that the b — > u amplitude prefers a small 

invariant mass of K + and K and hence a large invariant 
mass of K + and K~. In contrast, the b —> c amplitude 
prefers a small S23- The interference results to a minor 
hump shown in Fig.[^a). 



III. DECAY RATES AND CP ASYMMETRIES 



1 In our previous work |17|| we have argued that the spectrum 
should have a peak at the large m K + K - end. This is because we 
have introduced an additional nonresonant contribution to the 
lu- parameter parametrized as lj^ 11 = k 2p ^ P2 and employed 

the B~ — » D°K°K~ data and applied isospin symmetry to the 
B — * KK matrix elements to determine the unknown parameter 
k. Since this nonresonant term favors a small it^k+k s region, a 
peak of the spectrum at large m K + K - is thus expected. How- 



Results for the decay rates and CP asymmetries 
in B° — > K + K~ Kg( L ), K gKsKgr L \ are exhibited in 
Table [TIT] and Table [TV) respectively. We see that 



ever, such a bump is not seen experimentally [i| [see also Fig. 
[2jc)]. In this work we will no longer consider this term. 
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TABLE IV: Mixing-induced and direct CP asymmetries sin2/3 c ff (top) and Af (in %, bottom), respectively, in B° — > K + K~ Ks , 
KsKsKs and Kstv°iv decays. Results for (K + K~ Kl)cp± are identical to those for (K + K~ Ks)cp^- Experimental results 
are taken from H. 



Final state 


sin 2/3 eff 


Expt. 


(K K~ Ks)<t,K s excluded 

(K+K-K s )cp+ 

(K + K~ Kl)4,K l excluded 

KsKsKs 
KsKsK L 
K S TT°n 


n 791 +u.uuu+u.uui+u.u(jy 
u. /Zl_ ooi-0. 001-0. 020 

n 70C+0. 002+0. 007+0. 008 
u - ' ZD -0. 002-0. 004-0. 019 

n 701+0.000+0.001+0.009 

u - ,zl -0. 001-0. 001-0. 020 
II 71 O + 0. 000 + 0.000 + 0. 008 

u. 1 ia_ .000-0.000-0.019 

n 71 O+0. 000+0. 000+0. 008 
U. ( lO_ .000-0.000-0.019 
n 79Q+0. 001+0. 001+0. 009 
U. 1 ZJ_ 001-0. 001-0. 019 


0.58±0.13i^; 9 

0.58 ± O.lSto ll 
0.51 ± 0.21 

-0.84 ±0.71 ±0.08 




A f (%) 


Expt. 



(K K t\S)<pK s excluded 

(K+K-K s )cp+ 
{K + K-K L ) 

4>Kl excluded 

KsKsKs 
KsKsK L 
K s n°7v° 



-4. 
-4. 
-4. 

0. 

0. 

0. 



1 1 -t-l.Zzi-t-U.^l-t-U.^D 

ii -0. 91-0. 42-0. 30 
O7+1.30+0.45+0.37 
°'-l. 00-0. 46-0. 31 
1 1 +1.22+0.41+0.35 
11 -0. 91-0. 42-0. 30 
en+0. 01+0. 01+0. 05 
uy -0. 01-0. 01-0. 06 

Qn +o. 00+0. 00+0. 06 

yu -0. 00-0. 02-0. 08 

09+0.05+0.02+0.02 

oz -0. 03-0. 01-0. 03 



-15 ±9 

-15 ±9 
23 ± 15 

-27 ± 52 ± 13 



the predicted branching ratios are consistent with the 
data within the theoretical and experimental errors. 
From Table IIVI we can calculate the deviation of the 
mixing-induced CP asymmetry in B° — ► K + K~ Kg and 
KsKsKs from that measured in B — > (f> cE Ks, namely, 
Asin2/3 e ff = sm2(3 c g — sin 2fi^ c -K s ■ Using the fitted 
CKM's sin 2/3 = 0.695±g;gi| we obtain 

Asin2/3 if + K -K s = 0.026l°;°3o- 
Asin2/3 KsKsKs = 0.024±o ; ° 2 2 °, 
Asm2fe s7r07r o = 0.034±°;^o. (7) 



It should be stressed that despite the presence of color- 
allowed tree contributions in P>° — * K + K~ KgiL)> the 
deviation of the mixing-induced CP asymmetry in this 
penguin-dominated mode from that measured in B — > 
4> C cKs is very similar to that of the KsKsKs mode. 
However, direct CP asymmetry of the former, being of 
order 4%, is more prominent than the latter. 

Acknowledgments 

I'm grateful to Chun-Khiang Chua and Amarjit Soni 
for fruitful collaboration. 



[1] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 

71, 091102 (2005). 
[2] K. Abe et al. (Belle Collaboration), hep-ex/0507037 
[3] B. Aubert et al. (BaBar Collaboration) , |hep-ex/0507016| 
[4] See web page of Heavy Flavor Averaging Group, 

I htt p : / / www, slac . Stanford . edu/ xorg/hf ag | 
[5] D. London and A. Soni, Phys. Lett. B 407, 61 (1997). 
[6] Y. Grossman, Z. Ligeti, Y. Nir, and H. Quinn, Phys. Rev. 

D 68, 015004 (2003). 
[7] G. Engelhard, Y. Nir, and G. Raz, Phys. Rev. D 72, 

075013 (2005). 

[8] M. Gronau and J.L. Rosner, Phys. Lett. B 564, 90 
(2003). 

[9] B. Aubert et al. (BaBar Collaboration), |hep-ex/0607112| 
[10] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 

74, 032003 (2006). 
[11] A. Garmash et al. (Belle Collaboration), Phys. Rev. D 



71, 092003 (2005). 
[12] C. L. Y. Lee, M. Lu, and M. B. Wise, Phys. Rev. D 46, 
5040 (1992). 

[13] H. Y. Cheng and K. C. Yang, Phys. Rev. D 66, 054015 
(2002). 

[14] H.Y. Cheng, C.K. Chua, and A. Soni, in preparation. 
[15] S.J. Brodsky and G.R. Farrar, Phys. Rev. D 11, 1309 
(1975). 

[16] A. Garmash et al. (Belle Collaboration), Phys. Rev. Lett. 

96, 251803 (2006); Phys. Rev. D 75, 012006 (2007). 
[17] H.Y. Cheng, C.K. Chua, and A. Soni, Phys. Rev. D 72, 

094003 (2005). 

[18] CKMfitter Group, J. Charles et al, Eur. Phys. 
J. C 41, 1 (2005) and updated results from 
http://ckmfitter.in2p3.fr; UTfit Collaboration, M. Bona 
et ffl/., jhep-ph/050lT99"P 



